Crossover phenomena involving the dense O(n) phase 
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We explore the properties of the low-temperature phase of the O(n) loop model in two dimensions 
by means of transfer-matrix calculations and finite-size scaling. We determine the stability of this 
phase with respect to several kinds of perturbations, including cubic anisotropy, attraction between 
loop segments, double bonds and crossing bonds. In line with Coulomb gas predictions, cubic 
anisotropy and crossing bonds are found to be relevant and introduce crossover to different types 
of behavior. Whereas perturbations in the form of loop-loop attractions and double bonds are 
irrelevant, sufficiently strong perturbations of these types induce a phase transition of the Ising 
type, at least in the cases investigated. This Ising transition leaves the underlying universal low- 
temperature O(n) behavior unaffected. 
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I. INTRODUCTION 

The O(n) spin model is defined in terms of n- 
component spins on a lattice, with spin-spin interactions 
that satisfy O(n) symmetry, i.e., the model is isotropic 
in the n-dimensional spin- vector space. The cases n = 1, 
2 and 3 correspond with the Ising, XY and Heisenberg 
models respectively, but the significance of the O(n) 
model goes beyond these spin models. A loop expan- 
sion of the partition function of certain two-dimensional 
O(n) spin models 0, [|J leads to a system of noninter- 
secting loops, while the spin degrees of freedom are in- 
tegrated out. The resulting loop gas is called the O(n) 
loop model, sometimes abbreviated to just O(n) model. 
It has only discrete degrees of freedom, but the spin di- 
mensionality n appears in the partition sum of the O(n) 
loop model as a continuously variable parameter. In the 
limit n — > 0, the model serves to describe the behavior of 
polymer configurations 

For some two-dimensional O(n) spin models, a map- 
ping on a loop model is possible such that it yields the 
partition function in a form that enables the derivation of 
exact results [oT-fl3l|. These results show that there exist 
several "branches" of universality classes that continu- 
ously depend on the parameter n for — 2 < n < 2. 

One of these branches describes the phase transition 
between the high-temperature disordered spin phase and 
the low-temperature phase, where the spins display long- 
ranged correlations. In the terms of the loop model, the 
high-temperature phase is characterized by small loops 
and a low loop density, and the low-temperature phase 
by a high loop density and the existence of a loop of 
divergent size. 

The low-temperature phase appears to be more inter- 
esting than what one might expect on the basis of the 
known properties of the long-range ordered 0(1) or Ising 
model. For general n in the interval — 2 < n < 2, the 
low-temperature phase is still critical in the sense that the 
correlation functions display power-law behavior. Its uni- 
versal properties are described by another exactly solved 



branch. Recently, an exact transformation was applied to 
map the low-temperature branch of the O(n) loop model 
onto a tricritical loop model that includes vacant sites 
[Til . This mapping was applied for the case of the hon- 
eycomb as well as for that of the square lattice. Since 
this tricritical 0(n) model should have two more rel- 
evant temperature-like fields than the low-temperature 
branch, one may wonder whether these relevant direc- 
tions have some physical meaning in the low-temperature 
O(n) phase. 

The present work focuses on the stability properties of 
the low-temperature phase of the loop model with respect 
to several perturbations that move the loop model away 
from the exactly solvable point. These perturbations are: 

1. An attractive potential associated with loop seg- 
ments that collide at a vertex of the lattice; 

2. The introduction of double bonds, which allow 
some lattice edges to be covered by up to two loop 
segments; 

3. A cubic perturbation of the O(n) spin symmetry, 
which translates into the connection of four incom- 
ing loop segments at a vertex; 

4. Crossing bonds coupling O(n) spins, which cor- 
respond with crossing loop segments in the loop 
model, without affecting the O(n) symmetry of the 
corresponding spin model. 

The existing results in the literature, in particular from 
Coulomb gas theory [ll| HU predict, or at least suggest, 
the effects of these perturbations. Cubic deviations from 
O(n) symmetry were concluded to be irrelevant on the 
critical branch for n < 2, and to be relevant on the n < 2 
low-temperature branch [l6j]. Crossing bonds are pre- 
dicted to be described by the same exponent, so that they 
should also be relevant in the low-temperature phase. At- 
tractions between loop segments were however concluded 
[TtJ to be irrelevant in this phase. 
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Our present work purports to test the theoretical pre- 
dictions numerically, by means of transfer-matrix calcula- 
tions. In Sec.[TI]we define the models under investigation, 
and summarize the relevant existing results. Section IIIII 
explains the numerical procedures, and Sec. II VI presents 
the numerical results, concerning the phase diagram and 
the relevance or irrelevance of the various perturbations. 
We conclude with a discussion of the results in Sec. [V] 



II. MODELS 

The O(n) spin model with pair interactions is de- 
scribed by the reduced Hamiltonian 



u 



— y 



J(Si ■ Sj) , 



(1) 



where the sum is over all nearest-neighbor pairs, and the 
Si are n-dimensional vectors whose label indicates the 
site number i. They are normalized and 
their integration measure is J dst = 1. The function J 
describes the pair energy as a function of the spin prod- 
uct, and is usually chosen as a multiplicative constant, 
although other choices still preserve the O(n) symmetry. 
For the special choice J(y) = fceThi^l + zy) the Hamil- 
tonian becomes 



% = - ^ ln(l + zsi ■ sj) , 



(2) 



<ij> 



where the parameter z represents the coupling strength 
between neighboring 0(n) spins, and can thus be un- 
derstood as a measure of the inverse temperature. We 
consider the ferromagnetic case z > 0. The partition 
integral of this model can be written as 



Z = 



Y\_ J ds k 



Y[ (1 + ZSi ■ Si) 



(3) 



analysis appears possible for special values z = 

given by 



1/^2 ± v / 2~~^, 



- 2 < n < 2. 



(5) 



where the plus sign corresponds with a critical ordering 
transition separating the high-tcmperature phase from 
the low-temperature phase. The minus sign corresponds 
with the low-temperature O(n) phase. The solutions 
with the plus-sign are called branch 1, those with the 
minus sign branch 2. The exact results include the lead- 
ing scaling dimensions. 

These theoretical analyses are possible because of the 
special form of the spin-spin interaction J in Eq. ([TJ and 
because each spin occurs at most to the third power in 
the expansion of Eq. ©. 

For n > 2 the model described by Eq. Q no longer dis- 
plays critical points resembling branch 1 or 2, but there 
exists a line of critical points [18[ resembling the hard- 
hexagon transition [19| . 



B. Loop model on the square lattice 

Analogous to the case of the honeycomb lattice, an 
O(n) spin model can be defined such that it can be trans- 
formed into a system of nonintersecting loops [Til [l7| on 
the square lattice. In this case, the spins are located on 
the middle of the edges connecting the vertices of the 
0(n) loop model. The partition function of the latter 
model is a function of the loop weight n and the vertex 
weights u, v and w. The vertex weights of the square 
lattice O(n) model are defined in Fig. [T] 

The partition sum of the resulting loop model is simply 
written in terms of these weights as 



J <ij> 



Z\oop — y 



N„ N„ N,„ N, 

u v w n 



(6) 



where the products are on the sites and on the nearest- 
neighbor pairs of the lattice respectively. 



A. Loop model on the honeycomb lattice 

For the model on the honeycomb lattice, a graph ex- 
pansion Q of Eq. © expresses the partition function in 
terms of a sum over all configurations Q of nonintersect- 
ing loops on the edges of the honeycomb lattice: 



The sum is on all graphs Q consisting of nonintersecting 
loops on the square lattice, and N u , N v and N w are the 
numbers of vertices with weights u, v and w respectively. 

The resulting square-lattice O(n) loop model is solv- 
able for special choices of the vertex weights [TT| | . The 
solution includes four branches of critical points, where 
"critical" refers to algebraic decay of correlations. These 
four branches form a one-parameter family, parametrized 
by an angle 9 that is a four-valued function of the loop 
weight n. For branch k (with 1 < k < 4) the relation is 



J \oop 



* Nb n Nl 



W n = -2 cos(20) with ^> < 6 < , (7) 



where the graph Q covers Nb bonds of the lattice, and 
consists of Ni closed, nonintersecting loops. Each lattice 
edge may be covered by at most one loop segment. Exact 



with 



k' = 2, 1, 3, 4 for k = 1, 2, 3, 4 respectively. (8) 
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FIG. 1: The vertex weights u, v and w of the O(n) loop 
model on the square lattice. They are normalized such that 
the empty vertex has weight 1. In the models described by 
branches 1-4, the same weights apply to rotated versions of 
the vertices shown here. Also shown are additional vertices 
associated with perturbations of the model of Eq. ©, namely 
cubic vertices with weight c and crossing bond vertices with 
weight x. 

The vertex weights are 

4sin(0/2)cos(7r/4-0/4) 
U ~ ± 2 - {1 - 2sin(0/2)}{l + 2sin(<9/2)} 2 

„ = ± l + 2sin(fl/2) (9) 

2-{l-2sin(6»/2)}{l + 2sin(<9/2)} 2 V ; 
1 

w = . 

2 - {1 - 2 sin(0/2)}{l + 2 sin(6>/2)} 2 

It appeared that, after the relabeling of k! by k, branches 

1 and 2 share the universal properties of branches 1 and 

2 respectively on the honeycomb lattice. Branch 3 repre- 
sents a multicritical point where the O(n) critical tran- 
sition, a first-order transition, and an Ising transition 
merge. The Ising degrees of freedom can be understood 
in terms of dual spins on the faces of the square lattice, 
such that neighboring dual spins have the same sign only 
if they are separated by a loop segment (t7| . The uni- 
versal properties of branch 4 indicate a superposition of 
an Ising-like critical state and the low-temperature O(n) 
phase |17j | . Thus branch 4 is interpreted as a point where 
the aforementioned Ising degrees of freedom undergo an 
ordering transition. A sketch of the resulting phase dia- 
gram, as conjectured in Ref. and confirmed in Ref. ufj 
for the case n = 0, is reproduced in Fig. [2j 

The introduction of a sufficiently strong attractive po- 
tential between loop segments associated with the weight 
w, can, in principle, lead to an O(n) tricritical point [2lj. 
The latter result applies only to the case n — 0. This 
tricritical point is however of a different universal type 
as the multicritical point in the phase diagram of Fig. [5] 
Branches of tricritical points, parametrized by n, have 
been found for the square [lj, [HJ and the honeycomb 
[231 lattices with vacancies. 
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FIG. 2: Qualitative phase diagram of the O(n) loop model 
on the square lattice in a plane parametrized by the weights 
u ~ v and w. The locations of the exactly solved points 
are indicated by the symbol O and the corresponding branch 
number. The dilute O(n) loop phase corresponds with the dis- 
ordered phase of the O(n) spin model, and the dense phase 
with the low-temperature phase of the spin model. The full 
curve represents the critical line of the O(n) ordering tran- 
sition. The dense phase and the even denser Ising-ordered 
phase are separated by a line of Ising-like critical points, which 
is shown as a dashed line. The first-order transition between 
the disordered phase and the Ising-ordered low-temperature 
phase is shown as a double line. 

C. The n-component cubic model 

In the face-cubic spin models, the spin vector is re- 
stricted to lie along one of n Cartesian axes. Since it can 
still point in both directions of each axis, it has 2n pos- 
sible states. The spins lie on a lattice and have nearest- 
neighbor couplings of the form 

H = - Y^i Ks "i ■ s l + M ( s "i ■ S D 2 ] • ( 10 ) 

(ij) 

A graph expansion of the n-component cubic model was 
described in Refs. and O The resulting partition 
sum then depends, just as in the case of the 0(n) model, 
continuously on n. The form of the interaction be- 
tween the cubic spins is, as in Ref. [2|| chosen such that 
e M cosh K — 1 , in which case the graph expansion con- 
tains only even vertices, i.e., vertices that connect to an 
even number of neighbor sites. Thus, the graph expan- 
sion of the cubic model introduces a new vertex with 
four connected legs in comparison with the nonintersect- 
ing 0(n) loop model. The graph representation of the 
partition function of this cubic model is 

Z cnb = (2nfJ2^ Nl > (11) 

where z cu b = n~ 1 e AI sinh K, and N is the number of 
sites of the lattice, and the sum is over all graphs Q c that 
contain only even vertices, i.e., vertices connecting to 0, 2, 
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or 4 loop segments. The partition sum can be written in 
a form similar to Eq. ©, with u = v = z cu b, w = 0, and 
an additional cubic four-leg vertex with weight c = z ^ h . 

Also the Coulomb gas analysis of the O(n) model [16| 
uses four-leg vertices to describe a cubic perturbation. 
It predicts that cubic perturbations are irrelevant on the 
O(n) critical line for n < 2. Numerical analyses of the 
scaling dimensions of the n-component cubic critical line 
[24l [25j do indeed confirm O(n) universal behavior. In 
contrast, the cubic perturbation is predicted to be rel- 
evant in the low-temperature phase. It may thus seem 
rather curious that the cubic model of Eq. (fTTj) . when 
defined on the honeycomb lattice, reduces exactly to the 
form of Eq. (U), so that the cubic perturbation plays no 
role. The low-temperature phase of the n-component 
cubic model of Eq. (fTTj) on the honeycomb lattice is 
still in the universality class of the branch-2 O(n) model 
for n < 2. The exact results for the honeycomb O(n) 
model, including the critical point, apply as well to the 
n-component honeycomb cubic model of Eq. For 
the square-lattice cubic model the critical point is not 
exactly known. 



where the plus sign applies to branch 1 and the minus 
sign to branch 2. 

For the magnetic dimension Xh, one has mi = —7712 = 
1/2 and e\ = e2 = 1 — g, so that 



X,, 



2g 8 ' 



(14) 



In contrast with the iw-type vertex, which is believed to 
be irrelevant in the dense phase of the O(n) loop model 
fl7| , the cubic and the crossing-bond vertices change the 
topology of the graph representation. In the language of 
the mapping on the Coulomb gas, they are described by 
the four-leg watermelon diagram, which translates into 
magnetic charges mi — —m-i — 2 [5j. The scaling dimen- 
sions X c of a cubic perturbation of the O(n) symmetry, 
and X x of crossing bonds, are thus 



X r — X T — 1 — — -)- — . 

2g 2 



(15) 



For n = 2 these perturbations are marginal, and for n < 2 
they are relevant on branch 2. 

The temperature dimension is [l6j 



D. The 0(n) model with crossing bonds 

The perturbation of the low-temperature O(n) phase 
by the introduction of a square-lattice vertex with cross- 
ing bonds is predicted to be relevant @, [l6{ and is thus 
expected to introduce crossover to different universal be- 
havior. According to Jacobsen et al. [28|], the generic 
O(n) low-temperature phase is described by the crossing- 
bond model of Ref . [26|, [27]. This equivalence indicates the 
existence of a magnetic dimension Xh = 0, but attempts 
to verify this by finite-size scaling (2(| [28[ suffer from 
poor convergence, which may be attributed to logarith- 
mic factors. 



E. Exact exponents 

The exact results for the critical exponents of the crit- 
ical and the dense phase of the O(n) model can be con- 
veniently expressed in terms of the Coulomb gas scaling 
dimensions X(e.\, e2, mi, m^) associated with two pairs 
(ei,mi) and (e2,m2) of electric and magnetic charges. 
The scaling dimension associated with these two pairs is 

m 

A(ei,e 2 ,mi,m 2 ) = , (12) 

2g 2 

where g is the coupling constant of the Coulomb gas. For 
branches 1 and 2 of the O(n) model it is related to the 
loop weight n as 



X t = - 
9 



(16) 



which is irrelevant for n < 2 on branch 2. It is expected 
to describe the effects of a variation of the vertex weights 
u, v and w with respect to the branch-2 point. 



III. TRANSFER-MATRIX METHOD 

The transfer-matrix technique is used to calculate the 
free energy density and the magnetic correlation length 
of O(n) models wrapped on a cylinder of a finite circum- 
ference of L lattice units and of an infinite length. The 
free energy density of the system is 



f(L) = In A (L), 



(17) 



1 n 
= 1 ± — arccos— , 

7T 2 



(13) 



where Aq(L) is the leading eigenvalue of the transfer ma- 
trix T, and £ is the geometric factor, defined as the ratio 
of the unit of the finite size L over the layer thickness 
corresponding to the action of T. Thus, for the square 
lattice C = 1. 

It is useful to divide the transfer matrix into two diag- 
onal blocks or "sectors" as follows. When one cuts the 
cylinder through L edges which arc parallel to the axis 
of the cylinder, the number of dangling loop segments 
may be even or odd. It is obvious that the properties 
of evenness and oddness are conserved along the cylin- 
der, so that the transfer matrix decomposes into an odd 
and an even sector. The eigenvalue Aq(L) is the largest 
one in the even sector. The largest eigenvalue in the odd 
sector is denoted Ai(L). The states of the odd sector 
describe the effect of an additional single loop segment 
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running in the length direction of the cylinder. The map- 
ping between the O(n) spin model and the loop model 
provides the interpretation that the odd sector describes 
the spin-spin correlation function along the cylinder. The 
magnetic correlation length £h(L) is thus inversely pro- 
portional to the logarithm of the gap in the eigenvalue 
spectrum of T: 



^ 1 (£) = CMA (i)/Ai(L)]. 



(18) 



The calculation of the eigenvalues Aq(L) and Ai(£) of 
T is still made subject to the condition that the asso- 
ciated eigenvectors possess translational symmetry, i.e., 
the eigenvectors are invariant under the rotation of the 
cylinder over an angle lixjL about its axis. The trans- 
lational symmetry is in line with the form of the parti- 
tion sums given above and the periodic boundary condi- 
tions of a model on a cylinder. The correlation length 
^h(L) can be calculated numerically by the transfer- 
matrix method as a function of a parameter P repre- 
senting the distance (in some direction that remains to 
be specified) to a critical point or fixed point. Includ- 
ing this parameter in our notation, we define the scaled 
magnetic gap as 



X h {P,L) 



L 



(19) 



For models attracted by a conformally invariant fixed 
point, the scaled gap converges to the magnetic scaling 
dimension [30]. At a distance P of the fixed point, finite- 
size scaling [31| then predicts 



X h (P,L) = X h + aPL 



2-Xf 



(20) 



where X p is the smallest scaling dimension of the scaling 
fields to which P contributes, a is an unknown amplitude, 
and the dots stand for corrections to the leading scaling 
behavior that vanish for L — > oo. Since differentiation of 
Eq. ([201) yields 



dX h (P,L) 
dP 



aL 



2-Xf 



(21) 



it is possible to estimate Xp if numerical data for 
dXh(P, L)/dP are available for a range of finite sizes L. 
These data can be obtained by numerical differentiation, 
i.e., calculation of the scaled gap for several values of P, 
and subsequent fitting of a polynomial in P through the 
scaled gaps. 

The transfer-matrix construction for the O(n) model 
on the square lattice is described in Ref . [lj], including the 
coding that defines the transfer-matrix index in terms of 
the "connectivities" describing the topology of the loop 
configuration at a cross-section of the cylinder. A sparse- 
matrix decomposition allows the evaluation of the lead- 
ing eigenvalues of transfer matrices with linear sizes up 
to a few times 10 7 with the use of modest computer re- 
sources. Calculations for the O(n) model on the honey- 
comb lattice require a different sparse-matrix decompo- 
sition, which is explained in Ref. l29l 



Furthermore, we shall also introduce two new types of 
vertices on the square lattice that generate a larger set 
of connectivities than those of the nonintersecting loop 
model. These are the cubic vertex and the crossing-bond 
vertex, included in Fig. [TJ The introduction of the cu- 
bic vertex into the 0(n) loop model leads to connections 
between the loops and thus leads to a larger set of connec- 
tivities in comparison with the loop model. The number 
of mutually connected dangling edges is no longer re- 
stricted to two but may also assume multiples of 2. The 
coding and decoding needed for the construction of the 
transfer matrix corresponding with the cubic model of 
Eq. (TTT1) was described in Ref. H3. 

Also the presence of crossing-bonds leads to an in- 
crease of the number of connectivities. While the dan- 
gling loop segments can only be connected pairwise, the 
"well-nestedness" property of the non-intersecting loop 
model is lost. This property implies that, if dangling 
bonds i and j are connected, and dangling bonds k and 
I are also connected, that the situation i<k<j<l is 
excluded. Once crossing-bond vertices are allowed, the 
loops get entangled, and the situation i < k < j < I be- 
comes possible. The coding of this larger set of connec- 
tivities by means of integers 1,2, • • • is actually simpler 
than that of the well- nested connectivities [TtJ • The cod- 
ing is determined by a set of rules specifying an ordering 
of these connectivities. For completeness, we describe 
an ordering including the non-well-nested i-point con- 
nectivities. It is useful to represent a connectivity a by 
an array of integers i a = (i a (l), i Q (2), . . . , i a (L)), such 
that i a (k) = i a (l) if and only if the positions k and I are 
connected, i.e., if dangling edges k and I are covered by 
dangling segments of the same loop. The special value 
i m = represents a dangling edge not covered by a loop 
segment. The ordering of the connectivities, denoted by 
Greek symbols, is formulated in terms of these arrays of 
integers. The rules are: 

1. For the L Q -point connectivity a, remove the inte- 
gers with i m = from the array i a . This leads 
to an Lg-point dense connectivity without va- 
cancies, represented by an array Then, connec- 
tivity a precedes connectivity /? if exceeds the 
corresponding number Ls of connectivity /?. This 
provides only a partial ordering; it remains to or- 
der the set of connectivities with = L^. This 
remaining ordering will depend on the positions of 
the zeroes in i a , and on the connectivity a of the 
remaining dense configuration. 

2. Form an L„-bit binary number B a with (1) on 
position k if i a (k) ^ (i a (k) = 0). We can now 
specify that connectivity a precedes connectivity /3 
if L a — Lp and B a < Bp. 

3. The remaining task is to order the dense connec- 
tivities a. The first part is to find the position 
ria of the loop segment connecting to the loop seg- 
ment on position 1, i.e., the number that satisfies 
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*q(1) = ia( n &)- Then, we specify that dense con- 
nectivity a precedes j3 if n a < . 

4. If n a = rig, we define an L a > = L a — 2-point dense 
connectivity a' by removing positions 1 and n a 
from a. The remaining ordering is provided by the 
recursive application of the last two steps, adding 
primes at each new iteration, until a decision is 
found. 

The enumeration on the basis of this ordering requires 
only some trivial bookkeeping, involving numbers of con- 
nectivities of the relevant typ es, using methods presented 
already in Refs. [lj] and l32t An inverse algorithm that 
derives an array i a for a given connectivity number or 
transfer matrix index, was constructed similarly. 



IV. CROSSOVER AND THE DENSE PHASE 

We investigate the influence of various perturbations 
with respect to the branch-2 models on the honeycomb 
and the square lattice. 




-0.4 -0.2 0.2 0.4 0.6 0.8 1 



A. Attractions between loop segments on the 
square lattice 

We choose the plus signs in Eq. ([§]), while noting that 
the sign of u is irrelevant because the number of it-type 
vertices is even in the systems of interest. Since the num- 
ber of w-type vertices in a loop wrapping a cylinder with 
odd L is also odd, we have to keep in mind that the sign 
of v matters for odd system sizes. 

Denoting the vertex weights at the branch-2 and the 
branch-4 points as (u^ 2 \ v^ 2 \ w^) and (u^ 4 \ w^) 
respectively, we interpolate between the branch-2 and 
branch-4 points, and also extrapolate, by varying p in 



u(p) = (1 -p)u {2) +pu {A) 
v(p) = (l-p)v {2) +pv {i) 
w(p) = (1 — p)w^ +pw^ 



(22) 



We varied p in the range —0.5 < p < 1.5 and cal- 
culated the scaled magnetic gap Xh(p,L). Plots of this 
quantity as a function of p, for several values of the finite- 
size parameter L, are shown in Fig. [3] The main effect 
of increasing p is that the weight vo, which controls the 
attraction between neighboring loops, also increases. 

As a result, the loop configuration becomes denser, and 
at p — 1 it condenses into a state with Ising order as men- 
tioned in Sec. Ill Bl One observes that, for most n, there 
are clearly two different intersections of the curves for 
different L, near p = and p = 1. For p = the steep- 
est curves are those with the smallest L of the curves for 
different L. For p = 1 this situation is just the reverse. 
This shows that the perturbation with respect to branch 
2 is irrelevant, in agreement with the expected behavior 



FIG. 3: Scaled magnetic gap Xh versus the parameter p that 
interpolates between branch 2 (p — 0) and branch 4 (p — 1) 
for the O(n) model, (a) n = 0, (b) n = 0.5, (c) n = 1, (d) 
n = 1.5, and (e) n = 2. Results are shown for finite sizes 
L = 2, 4, • • • , 12. The slope of the curves increases with L on 
the right-hand side. 



TABLE I: Numerical estimates of the scaling dimensions X v 
and X W2 , which belong to the most relevant scaling fields to 
which p and W2 contribute respectively. 



n 


Xp 







1.000 (1) 


1.0(1) 


0.5 


1.0000 (1) 


0.99(1) 


1 


1.000000 (2) 


0.9999(1) 


1.5 


1.00001 (1) 


1.001(1) 


2 


1.01 (1) 


0.999(1) 



for the leading thermal exponent according to Eq. (|16[) . 
In contrast, the perturbation due to the variation of p 
with respect to the branch-4 point is seen to be relevant. 
Numerical differentiation to p of the scaled gaps at the 
branch-4 points for several n, and subsequent analysis 
according to Eq. (l2~Tj) (with p instead of P) yielded esti- 
mates of X p that are shown in Table HI One observes that 
the results for X p are close to the known Ising tempera- 
ture dimension Xt = 1. This confirms the Ising nature of 
the transition driven by p, in line with the conclusion 17 1 
that it takes place independent of the critical background 
of the dense O(n) model. 
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B. The nature of the Ising-ordered phase 

Figure [3] in the preceding subsection, and numerical 
results for larger values of p > 1, indicate that the scaled 
gaps increase approximately linearly with the finite size, 
and thus that the magnetic correlation length becomes 
constant. This corresponds with a magnetic correlation 
function that decays exponentially in the infinite plane. 
In this respect, the dense phase and the Ising ordered 
phase, separated by the Ising line as shown in Fig. [2J 
are different. Still, these two phases are assumed [171 ] 
to share the basic universal properties of the dense O(n) 
loop model. We test this assumption by eliminating the 
reason why the scaled gaps increase sharply at p = 1 
and beyond. The reason is that, near the line of Ising 
transitions, the loop configurations become so dense that 
most vertices are of the w-type. Since the even sector of 
the transfer matrix allows only loop configurations that 
cover an even number of edges in the transfer direction, 
only configurations of the even sector fit well on a lattice 
with even L. This explains the increase of the gap be- 
tween the even and the odd sector as the the Ising line is 
approached. Similarly, such dense configurations in the 
odd sector will only fit well on lattices with odd L. 

In order to define a type of magnetic gap that excludes 
these effects of even-odd alternation, one has to select 
even or odd systems in accordance with the sector. Wc 
thus define a magnetic scaled gap X m (p,L) as 



X m (p,L) 
for odd L, and 



CL VAo(X + l)Ao(L-l) 

2^ ln ML) (23) 
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FIG. 4: Scaled magnetic gap X m {p,L) versus the parameter 
p that interpolates between branch 2 (p = 0) and branch 4 
(p — 1) for the cases (a) n = 0, (b) n = 0.5, (c) n = 1, 
(d) n = 1.5, and (e) n = 2 O(n) model. Results are shown 
for L = 3, 4, • • • , 12. The dashed curves with symbols x 
are used for even L, and the solid curves with symbols + for 
odd L. Both sets of curves tend to converge to a common 
p-independent limit when L increases. 



2tt y/K^L + 1)^-1) 



(24) 



for even L. 

The results for these scaled magnetic gaps are shown 
in Fig. |4] as a function of p, for several values of n. These 
results indicate that, for n < 2, the universal character of 
the magnetic correlations in the dense phase is indepen- 
dent of p. In particular, it remains unchanged under the 
Ising transition and the onset of the Ising-type long-range 
order. 

Next we estimated the associated magnetic scaling di- 
mension X m at the Ising point (branch 4) by fitting the 
numerical results using Eq. (|20[) for several n. The results 
are listed in Table |TTJ Here our choice of the sign of v in 
Eq. dU follows the change of sign of w such that vw > 0. 
The change of sign of w is caused by the change of sign of 
the common denominator in Eq. ([§]) , which arises because 
the weight of the empty vertex is normalized to +1. A 
consequence of this change of sign is that there will be a 
jump with value 1/8 in X m near n — 0.087378025, where 
the change of sign of w occurs, if the weight v is kept 
positive. The jump is equal to the "interface dimension" 
denoted X- lnt i in Ref. Uji 



TABLE II: Numerical results for the magnetic scaling dimen- 
sion X m of branch 4 (rightmost column), compared to the 
exact results for the magnetic dimension X% for branch 2 and 
for branch 4. The branch-4 magnetic dimension is equal to 
the branch-2 dimension plus 1/8. 



n 


Xh (branch 2) 


X h (branch 4) 







-0.1875 


-0.0625 


-0.18749 (1) 


0.075 


-0.16865534 


-0.04365534312 


-0.168653 (2) 


0.0875 


-0.16561812 


-0.04061812640 


-0.16562 (2) 


0.1 


-0.16260929 


-0.03760929546 


-0.162610 (2) 


0.125 


-0.15667507 


-0.03167507924 


-0.156675 (2) 


0.5 


-0.07909087 


0.04590912236 


-0.07909 (1) 


1 





0.125 


0.0000000 (1) 


1.5 


0.061874313 


0.18687431332 


0.0618743 (1) 


2 


0.125 


0.25 


0.125000 (1) 



C. Double bonds in the honeycomb model 

The loop model of Eq. Q is extended by allowing the 
edges with one of the three possible orientations, say the 
vertical edges, to be covered by up to two loop segments. 
The honeycomb lattice can be decomposed in building 
blocks consisting of a vertical edge and the pair of vertices 
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/--/ x 



z ~~ v 



-w 



~w 




FIG. 5: The vertex weights of the 0(n) loop model on the 
square lattice as obtained by a mapping on the honeycomb 
lattice. These vertices allow the occupation of vertical edges 
of the honeycomb lattice by two loop segments. The resulting 
vertex weights allow for a factor z per loop segment on the 
honeycomb lattice. 



at its ends. These units are shown in Fig.[5j together with 
their weights, which include a factor z per loop segment. 

As indicated in Fig. [51 each such unit of the honey- 
comb lattice can be replaced by a square-lattice vertex. 
This substitution maps the honeycomb model with dou- 
ble bonds on an O(n) loop model on the square lattice, 
but the vertex weights, also shown in Fig. [5j are not of 
the form of Eq. © because the weights of the u- and w- 
type vertices depend on their orientation. Expressed in 
the enlarged set of square-lattice weights, the partition 
sum takes the form 



7 _ ^"1 

■^■loop — / J u \ 

Q 



(25) 



where the indices appended to u and w indicate the orien- 
tation of the vertex. The vertex weight z' — W2 describes 
a double bond covering a lattice edge. The resulting con- 
nectivities are, however, still of the nonintersecting loop 
type, which means that the two loop segments on an edge 
do not cross or mutually connect. 

We evaluated the scaled gaps according to Eqs. (fT8|) 
and (|19p of the model of Eq. (f2"5"f for several values of 
the finite size L, using the weights u±, 1*2, V and w\ as 
obtained from the equivalence with the branch-2 point 
of the honeycomb model with z according to Eq. (J5j. 
The influence of double bond bonds was determined by 
including several nonzero values of the weight z' = W2- 
The calculation of the scaled gaps used the geometric 
factor £ = 2/V3 for the honeycomb lattice. Since the 
additional weight z' introduces anisotropy, the asymp- 
totic conformal symmetry is broken, and the scaled gaps 
for z' 7^ no longer directly relate to the scaling dimen- 
sion Xfr. For this reason we add a tilde and denote the 
scaled gaps as Xh(z', L). The results are shown in Fig.|6] 

While the interpretation of the scaled gap in terms 
of the scaling dimension Xh is no longer valid, the in- 
tersections on the right-hand side of these figures, with 
slopes increasing with L, still indicate that a phase tran- 
sition takes place, resembling the Ising-like ordering for 
the square lattice model in Sec. IIV Al The Ising charac- 
ter of this transition was verified by means of numerical 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

FIG. 6: Scaled magnetic gap Xh(w2, L) versus the weight W2 
of a double bond covering an edge of the honeycomb lattice 
for O(n) model with (a) n = 0, (b) n = 0.5, (c) n — 1, (d) 
n = 1.5, (e) n — 2. Results are shown for finite sizes L = 2, 
4, ■ • ■ , 12. The slope of the curves increases with L on the 
right-hand side. 



differentiation of Xh(w2, L) with respect to W2 in the in- 
tersection points and finite-size scaling, analogous to the 
analysis of X p in Table Q] The numerical estimates of the 
scaling dimension X W2 for several n are included in Ta- 
ble HI They are close to the Ising temperature dimension 
X t = l. 



D. Cubic anisotropy 

The introduction of cubic vertices into the 0(n) model 
of Eq. ([S]) modifies the partition function as follows 



J \oop 



EN„ AT,, N„, N r Ni 
u u v v w w c °n 1 



(26) 



where N c is the number of vertices of type c. For w = 0, 
v = u and c = u 2 it reduces, apart from a multiplica- 
tive constant, to the partition sum of the cubic model, 
Eq. (fTTj) . Thus, Eq. (|2"fJ)) can interpolate between the 
nonintersecting loop model and the cubic model. 

We investigate the effect of cubic perturbations by 
varying c, while keeping the other vertex weights fixed 
at their branch-2 values. The results for the scaled gaps 
are shown in Fig.[7J For n — 0, the vanishing loop weight 
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TABLE III: Scaling dimensions X c associated with a cubic 
perturbation, and X x associated with crossing bonds for the 
case of the dense O(n) phase. These results are obtained by 
numerical differentiation of the scaled magnetic gap Xh(L) 
at the branch 2 point. The vanishing of some numerical re- 
sults for n — and 1 is an artefact due to suppression of the 
respective critical amplitudes. 





Coulomb gas (branch 2) 


numerical results 


n 


9 


X x — X c 


x c 


X x 





0.5 


0.75 





0.75001 (1) 


0.5 


0.58043 


1.00922 


1.009 (1) 


1.0092 (2) 


1 


2/3 


1.25 








1.5 


0.76995 


1.50552 


1.505 (1) 


1.5055 (3) 


2 


1 


2 


2.000 (2) 


1.998 (2) 



gas predictions. 



0.1 0.2 0.3 0.4 0.5 



E. Crossing bonds 



FIG. 7: Scaled magnetic gap Xh versus the weight c of cubic 
vertices introduced into the dense phase of the O(n) noninter- 
secting loop model according to Eq. (|26[) . Results are shown 
for finite sizes L = 2, 4, • • • , 12. (a) n = 0. The weight c 
has no effect because the loop weight n — does not allow 
a nonzero density of c-type vertices. The gap increases as a 
function of L. The cases n = 0.5, 1, 1.5, 2 are shown in (b)-(e) 
respectively. The slope of the curves increases with L on the 
right-hand side. 



prevents the introduction of cubic vertices, so that the re- 
sults do not depend on c. The slopes of the curves with 
< n < 2 are seen to increase with L near c = 0, which 
shows that the cubic perturbation is relevant. Further- 
more, some of the plots display two more sets of intersec- 
tions, of which the middle ones are indicative of a stable 
fixed point, and the rightmost ones of an unstable fixed 
point, resembling that of the Ising-like transition induced 
by the w-type vertices. These results are consistent with 
the interpretation that, in the range attracted by the sta- 
ble fixed point, the scaled gap converges to Xh — for 
< n < 2 in a range c > 0. This could be confirmed by 
numerical extrapolations at some values of c, all of which 
satisfied X h < 0.01. 

In order to numerically determine the exponent re- 
sponsible for the cubic crossover, we performed a numer- 
ical differentiation of the scaled magnetic gap Xh(c,L) 
with respect to the weight c of the cubic vertex at the 
point c = for several values of n. The finite-size data for 
this derivative were subsequently analyzed according to 
Eq. (|2l"j) . with P replaced by c. The resulting estimates 
of X c are shown in Table IIII1 together with the Coulomb 



We next introduce, starting from the branch-2 low- 
temperature points of the square-lattice loop model, a 
nonzero weight of the crossing-bond vertex. We calcu- 
lated the scaled gaps for a range of values of the weight 
x, and for several values of n. The results are shown in 
Fig.El These data indicate that, for n < 2 and a range of 
x > 0, crossover occurs to a different universality class of 
dense intersecting loop models, with a magnetic exponent 
that is different from that of nonintersecting loop models 
for n ^ 1 . This interpretation is in line with a prediction 
of Jacobsen et al. [28| in terms of exact results (2(| [53] ■ 
We have attempted to find the conformal anomaly c a and 
the magnetic exponent Xh from the finite-size data for 
a few points in the phase diagram. Estimates for c a are 
obtained by fitting f(L) = /(oo) +7rc a /(6L 2 ) to transfer- 
matrix results for the free energy using three subsequent 
values of L. These estimates display slow apparent con- 
vergence and are thus hard to extrapolate. Extrapolation 
was done assuming finite-size dependence as L~ 2 . The 
results are listed in Table IIVI As a tentative error mar- 
gin we quote ten times the difference between the last 
two extrapolations. Also the data for Xh were hard to 
extrapolate; we simply quote the result Xh(L) obtained 
from Eq. (|19l) for L — 14, with a tentative error margin 
of 10 times X h (U) - X h (12). 

We also performed numerical differentiations of the 
scaled magnetic gaps Xh(x, L) with respect to the weight 
x of the crossing-bonds vertex at the point x = for sev- 
eral n. The finite-size data were subsequently analyzed 
according to Eq. (|2"Tj) . with P replaced by x. This yielded 
estimates of X x that are shown in Table lTnl together with 
the Coulomb gas predictions. 
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TABLE IV: Numerical estimates of the conformal anomaly c a 
and the magnetic exponent Xh in the dense O(n) phase with 
crossing bonds. 



n 


X 


Ca 


x h 





0.5 


-1.42 


(26) 


-0.066 


(29) 





0.8 


-1.36 


(29) 


-0.053 


(15) 


0.5 


0.4 


-0.64 


(19) 


-0.036 


(13) 


0.5 


0.6 


-0.65 


(10) 


-0.031 


(9) 


1 


0.8 


-0.01 


(10) 


0.000 


(2) 


1 


1.1 


0.05 


(64) 


0.011 


(50) 


1.5 


0.6 


0.55 


(2) 


0.031 


(3) 


1.5 


0.8 


0.53 


(11) 


0.025 


(3) 


1.5 


1.0 


0.52 


(35) 


0.025 


(30) 


2 


0.3 


1.000 


(2) 


0.095 


( 1) 


2 


0.6 


0.998 


(46) 


0.067 


( 1) 


2 


0.9 


0.97 


(29) 


0.046 


(35) 
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FIG. 8: Scaled magnetic gap Xh versus the weight of the 
crossing-bond vertex, for the cases (a) n = 0, (b) n — 0.5, (c) 
n — 1, (d) n — 1.5, (e) n = 2. Results are shown for finite 
sizes L = 2, 4, • • • , 12. The slope of the curves increases with 
L on the right-hand side. 



V. DISCUSSION 



The numerical results presented in Sec. IIV Al for the 
effect of the w-type vertex, representing loop-loop at- 
tractions, on the low-temperature phase of the square- 
lattice model, agree with the predicted [13] behavior. As 
indicated by the curves in Fig. |31 these attractions are 
irrelevant in the whole range interpolating between the 
branch-4 and branch-2 points. In contrast, they are rele- 



vant in the equivalent tricritical O(n) model with vacan- 
cies [l4j. The latter model is obtained by summing out 
part of the loops, which thus yields a system that is far 
more susceptible to attractions between the loops. 

Furthermore, the numerical results confirm that the 
loop-loop attractions are also relevant at the branch-4 
point, and that the transition that takes place at this 
point is Ising-like. Apart from the explanation of this 
transition in terms of the onset of long-range order of 
the dual Ising spins, it may be worthwhile to mention 
that the type of phase diagram in Fig. [3J including the 
Ising line, is reproduced by an Ising model with vacan- 
cies quoted in Ref. [IJ In the latter case, the Ising line 
corresponds with the onset of phase separation between 
a phase dominated by Ising spins and a phase dominated 
by vacancies. 

The analysis of the O(n) magnetic dimension X m pre- 
sented in Sec. IIV Bl confirms that the low-temperature 
0(n) universal character remains unaffected across the 
Ising transition on branch 4. This result was obtained by 
means of a careful formulation of the correlation function 
associated with X m , which takes into account the even- 
or oddness of the system, even in the limit of infinite size. 

On the basis of results [26l - [28j j for loop models that al- 
low crossing bonds, doubts have arisen to what extent the 
behavior found for the exactly solvable models of branch 
2 of the honeycomb Q , the square 11 1 , and the triangu- 
lar lattice [331 ] is representative for the low-temperature 
O(n) phase. Since multiple bonds arise in a natural way 
in graph expansions of more general spin O(n) models, 
we studied the effect of double bonds in Sec. IIV CI This 
perturbation was found to be irrelevant for branch 2 of 
the model on the honeycomb lattice. Its effect appears to 
be very similar to that of loop-loop attractions due to the 
type-w; vertex in the case of the square lattice. The sim- 
ilarity includes the Ising-like transition that takes place 
at a sufficient weight of the double bonds. 

The Coulomb gas prediction that the cubic perturba- 
tion is relevant for n < 2 is quantitatively confirmed by 
the results for X c in Table UTT1 only for n = 0.5 and 1.5. 
We attribute the vanishing numerical results for n = 
and 1 to vanishing amplitudes associated with the cubic 
perturbation. The numerical result X c = at n = is 
due to the fact that the zero loop weight excludes type-c 
vertices even at nonzero fugacity. For n = 1, the loop 
weights are equal to 1, so that the distinction between 
type-z and type-c vertices disappears, and so do the am- 
plitudes associated with the cubic perturbation. Taking 
into account this explanation for the deviating results 
for n — and 1, we conclude that there is a satisfactory 
agreement with the Coulomb gas predictions. 

As mentioned in Sec. Ill CI the relevance of cubic per- 
turbations in the low-temperature phase of the square- 
lattice model of Eq. (fTTj) for n < 2, which was confirmed 
in Sec. IIV Dl seems very peculiar in relation with the ab- 
sence of crossover to cubic behavior for the honeycomb 
model. In view of the mutually incompatible values of 
X] r , the conclusion that the dense phases of the partition 
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sum of Eq. pT|) display different universal behavior for 
the square and honeycomb lattices is inescapable, and 
has to be attributed to the low coordination number of 
the honeycomb lattice. 

One may wonder if a similar paradox occurs in the spin 
representation of the model, for which our physical intu- 
ition may provide further insight. The equivalence with 
Eq. (|11[) applies only to cases where n is a positive inte- 
ger. The case n = 1 fails to provide more clarity because 
the amplitude associated with c vanishes. For n = 2wc 
do expect a phase transition to a long-range-ordered state 
in the spin model described by Eq. (ITU)) at sufficiently low 
temperatures, but the condition e M cosh if = 1 excludes 
this low-temperature range from the loop representation 
of Eq. (|11[) . There is no apparent conflict with the ex- 
pected behavior of the cubic spin model. 

The introduction of the cubic vertices into the loop 
model on the square lattice in Sec. IIVDI yielded results 
that are consistent with the interpretation that the scaled 
gaps converge to for n < 2 in a range c > 0. This in- 
terpretation is in line with the expected long-range order 
of spin models with a cubic perturbation at low temper- 
atures. 

Furthermore we note that the mapping of the hon- 
eycomb model on the square lattice model presented in 
Sec. lIV Cl ean be extended to include cubic vertices on the 
square lattices, which corresponds to mutually connect- 
ing double bonds on the vertical edges of the honeycomb 
model. Thus, cubic crossover will occur on a suitably 
generalized honeycomb model. 

Also in the case of perturbations introduced by the 
crossing-bond vertex we find a satisfactory agreement 
with the Coulomb gas predictions. Also in this case the 
amplitude a in Eq. (|20l) due to the perturbation vanishes 
at n = 1 (see Table which explains the vanishing 
of the corresponding numerical result. For n — 1.5 the 



amplitude is still rather small (see Fig.(5fd)), but the nu- 
merical differentiation method is sufficiently sensitive to 
determine the scaling dimension X x . 

The scaled gaps in the dense phase perturbed by cross- 
ing bonds appear difficult to analyze. While our range of 
finite sizes is insufficient for reliable extrapolations, the 
data seem in line with Xh = and the occurrence of log- 
arithmic factors as predicted for crossing-bond models 
f26T - |28j ] . Moreover, the results for the conformal anomaly 
given in Table ITVl show a trend consistent with c a = n — 1 
as predicted by Martins et al. [2a |. 

In conclusion, our results confirm the phase diagram 
of the 0(n) model as conjectured in Ref. [l7|, in particu- 
lar the Ising transition between the dense and the Ising- 
ordered phases, and its Ising scaling dimension. Cubic 
anisotropy and crossing bonds are proved numerically to 
be relevant and introduce crossover to different universal 
behavior in the low-temperature (dense) phase. 

Finally, we note that our numerical results display so- 
called nonuniversal behavior as a function of the various 
perturbations in the n — 2 models, as expected from the 
mapping between the cubic model and, e.g., the Ashkin- 
Teller model and the equivalence of the latter model 
with the eight- vertex model [341-136} . 
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